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Abstract 

We study the mixing time of random graphs in the d-dimensional toric unit cube [0, l] d 
generated by the geographical threshold graph (GTG) model, a generalization of random ge- 
ometric graphs (RGG). In a GTG, nodes are distributed in a Euclidean space, and edges are 
assigned according to a threshold function involving the distance between nodes as well as 
randomly chosen node weights, drawn from some distribution. The connectivity threshold for 
GTGs is comparable to that of RGGs, essentially corresponding to a connectivity radius of 
r = (log n/n) 1 ^. However, the degree distributions at this threshold are quite different: in 
an RGG the degrees are essentially uniform, while RGGs have heterogeneous degrees that de- 
pend upon the weight distribution. Herein, we study the mixing times of random walks on 
d-dimensional GTGs near the connectivity threshold for d > 2. If the weight distribution func- 
tion decays with P[W > x] = 0(\/x d+u ) for an arbitrarily small constant v > then the 
mixing time of GTG is 0(n 2 / d (log n)( d ~ 2 )/ d ). This matches the known mixing bounds for the 
d-dimensional RGG. 

Keywords: Geographical threshold graphs, Mixing time, Random geometric graphs. 

1 Introduction 

In recent years, we have witnessed the development of numerous approaches to study the structure 
of large real-world technological and social networks, and to optimize processes on these networks. 
Large networks, such as the Internet, World Wide Web, phone call graphs, infections disease con- 
tacts and financial transactions, have provided new challenges for modeling and analysis [8J. As 
an example, Web graphs may have billions of nodes and edges, which implies that processing and 
extracting information on these large sets of data, is 'hard' pQ. Extensive theoretical and experi- 
mental research has been done in web-graph modeling, attempting to capture both the structure 
and dynamics of the web graph [HI [21 [231 [^1 ESI O El [TT] . 

In general, a particularly fertile approach is to consider the network as an instance of an ensem- 
ble, arising from a suitable random generative model. Since the seminal papers on the evolution of 
uniform random graph model [201 El], many other models have been proposed to better capture the 
structure seen in real- world networks, which are systematically covered in [19J. One straightforward 
example is the random geometric graph (RGG) model, where nodes are placed uniformly at ran- 
dom in a Euclidean space and edges are placed between any two nodes within a threshold distance. 
For further study of RGGs, see the monograph by Penrose [29J. The RGGs have the advantage 
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of describing many aspects of systems such as sensor networks, while avoiding unnecessary details. 
However, they fail to capture heterogeneity in the network. 

Geographical threshold graphs (GTGs), introduced in [28], are a generalization of RGGs. Het- 
erogeneity in the network is provided via a richer stochastic model that nevertheless preserves much 
of the simplicity of the RGG model. GTGs assign to nodes both a location and a weight. The 
weight may represent a quantity such as transmission power in a wireless network or influence in a 
social network. Edges are placed between two nodes if a symmetric function of their weights and 
the distance between them exceeds a certain threshold |13| . 

Structural properties of GTGs, such as connectivity, clustering coefficient, degree distribution, 
diameter, existence and absence of the giant component, chromatic number have been recently 
analyzed [11[ [T2"| 114] . These properties are not merely of theoretical importance, but also play 
an important role in applications. In communication networks, connectivity implies the ability to 
reach all parts of the network. In packet routing, diameter gives the minimal number of hops needed 
for transmission between two arbitrary nodes. In the case of epidemics, the existence or absence 
of the giant component controls whether the epidemic spreads or is contained. When treating the 
node colors as the different radio channels or frequencies, the chromatic number gives the minimal 
number of channels needed so that neighboring radios do not interfere with each other. 

Herein, we consider random walks on GTGs near the connectivity threshold. Random walks (or 
more formally, Markov chains) on large networks have many applications. For example, random 
walks model the spread of disease or the dispersion of information [10J. The mixing time of a 
random walk is the expected number of random steps that are required to guarantee that the 
current distribution is close to the stationary distribution. Mixing times are an essential tool in 
both theory and practice: for example, see the recent survey of Diaconis [18] on Markov chain 
Monte Carlo methods. In [6], the authors derived the mixing time of exponential random graphs, 
a model extensively used in sociology, showing that the mixing time of the Glauber dynamics was 
G(n 2 logn) for the unimodal Gibbs distribution, and exponential for the multimodal case. For the 
definitions of the Glauber dynamics and Gibbs distribution see [6]. 

Upper bounds on the mixing time for RGG at the connectivity threshold have been well studied. 
For the 2-dimensional RGG, Avin and Ercal [4] showed that the mixing time is 0(n). More 
recently, Cooper and Frieze [16] proved that for d > 3, the mixing time of a ^-dimensional RGG 
is (9(n 2//d ) (in this notation, the logarithmic factors are suppressed). In this paper, we study the 
mixing times of random walks on d-dimensional GTGs near the connectivity threshold, where 
d > 2. We prove that when the node weight distribution decays sufficiently quickly, the mixing 
time is 0{n 2 / d '(log n) (d_ 

-2)/ d ) = 0(n 2 l d ), which matches the mixing bounds for RGG. This result is 
formulated more precisely as Theorem [2] in the next section. 

2 The GTG model and the mixing time 

The GTG model is constructed from of a set of n nodes placed independently uniformly at random 
into the unit cube in ~R, d . A non-negative weight W{, taken randomly and independently from a 
continuous probability distribution function f(w) : R + — > R + , is assigned to each node Vi for 
i G {1,2,... , n}. For nodes i,j at distance r^, the edge {i,j} exists if and only if the following 
connectivity relation is satisfied: 

G{w u w 3 )D{r l3 )>9 n . (1) 

Here G(u)i,Wj) is the interaction strength between nodes, D{ r) is a decreasing function of r and 9 n 
is a given threshold parameter that depends on the size n of the network. The interaction strength 
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G(wi,Wj) is usually taken to be symmetric and either multiplicatively or additively separable, i.e., 
in the form of G(wi,Wj) = g(wi)g(wj) or G{wi,Wj) = g(wi) + g(wj). We use D (r) = r~ s where 
s > 0, which is a typical attenuation in the path-loss model in wireless communications [13] . 

Some basic results have already been shown, including the expected degree of a node with 
given weight w, when the nodes are distributed uniformly over a unit space [281 [13j. In both the 
multiplicative and additive cases of G(wi,Wj), questions of diameter, connectivity, and topology 
control have been addressed [13]. 

Here we restrict ourselves to nodes distributed uniformly over [0, l] d . For analytical simplicity 
we take the space to be the d-dimensional toric unit cube [0, l] d . Our connectivity relation uses an 
additive interaction strength G(w\,W2) = W\ + w>2 and a decay function D(r) = r d , so that nodes 
i,j are adjacent when 

> po 

This connectivity relation identifies a d-dimensional sphere of influence for each vertex. 

We assume that our weight distribution f(w) has finite mean and finite variance. Let the 
cumulative density function (cdf), for the distribution of node weights f(w), be 

F(x) = P[W < x] = [ f(w)dw. (3) 
Jo 

The argument in [12] characterizing the degrees of a GTG for a 2-dimensional GTG is easily 
generalized to dimension d (only the leading constant changes). For any threshold 6 n = 0(n) and 
any weight distribution such that P[W > x] = 0(l/x 1+e ) for an arbitrarily small constant e > 0, 
the degree distribution of a node v with weight w follows the binomial distribution 



&eg(v\w) ~ Bin (n — l,p(w)) , (4) 



where p(w) = -^(w + //), n = E[W] is the expected node weight and is the volume of the unit 
ball in d dimensions: 

TT d / 2 f ir k /k\, d = 2k even, 

d = T(d/2 + 1) = \ 2 d kl7r k /dl, d = 2k + l odd. 

Herein, we assume that F[W >x} = 0{l/x d+v ) = o(l/x 2 ) which ensures that the weight distribu- 
tion has finite mean and finite variance. 

We now highlight the differences between the GTG and RGG models. The main characteristic 
of the GTG model is its tunable topology. By changing the input parameters f(w) and n , one 
can obtain graphs with different structural properties. For example, we can generate an RGG for 
a desired degree distribution, while the degrees of RGG are always uniform. The major distinction 
between our analysis for GTGs and the analysis for RGGs, lies in addressing the following two 
issues: (i) two spatially close nodes in a GTG are not necessarily connected, since they may both 
have very low weight; (ii) two distant nodes in a GTG are not necessarily disconnected, since one 
of them can carry heavy weight. These two issues never happen in the RGG model, and represent 
a challenge in our analysis. 

Theorem 5.3 of |11] characterizes the connectivity threshold for a GTG in 2 dimensions. We list 
the changes to the proof to generalize to d-dimensional GTGs. We tile the unit space [0, l} d into 
0(n/logn) cubes of equal volume (as opposed to 2-dimensional squares). We use the connectivity 
relation equation ([2]) for the d-dimensional space and connectivity radius r conn = (log(cm) / (anT^)) 1 ^ 
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instead, with the constant a £ (0,1) specified in Theorem 5.3 of With these modifications, 
we obtain the proof of the connectivity threshold for a ci-dimensional GTG, given by the following 
theorem. 

Theorem 1. Let G be a GTG in the d- dimensional toric unit cube [0, l] d with threshold function 
9 n = cn/logn where the constant c < sup^g^i) ceF~ 1 (l — a)/A. Then G is connected wh]^. 

In essence, the proof of Theorem Q] consists of two parts. First, one shows that the an nodes of 
highest weight are connected. Next, one shows that the remaining (1 — a)n nodes are connected 
to first set. This partition into high weight nodes and low weight nodes will also be useful herein. 

This paper gives an upper bound on the mixing time for a simple random walk on a geograph- 
ical threshold graph G = (V(G),E(G)) at the connectivity threshold, provided that the weight 
distribution decays at an adequate rate. A simple random walk on a graph G (cf. [27]) consists of 
a sequence of vertices (u>n,ioi, . . . ,wt, . . ■) such that for t > 0, P[u;t + i = j | Wt = i] is 1/ deg(i) if 
{i,j} £ E(G) and otherwise. Let ir denote the stationary distribution of this random walk, so 
that 7Tfc = deg(k)/2\E\ for every node k. Supposing that % is the initial node of our random walk, 
let P t {i, •) denote the distribution of the states at time t. The variational distance at time t is 



A;(i) = ~£|P f (*,j) 



2 

jev 

When G is not bipartite, we have lim^oo Aj(i) = for every i 6 V{G). The mixing time from 
node i measures how quickly P (i, •) converges to n. Explicitly, this mixing time from i is defined 
as 

Ti{5) = min{t | A;(i') < 6,W > t}. 

The mixing time of G is 

t(5) = max tAS). 

i€V(G) 

We choose S = 1/n as our desired distance from the stationary distribution. Our main result is 
given by the following theorem. 

Theorem 2. Let G be a connected GTG with threshold function 6 n = cnj log n in the d-dimensional 
toric unit cube [0, l] d . If the weight distribution satisfies F[W > x] = 0(l/x d+u ) for v > 0, and 

c < - — -mini sup aF~ 1 (a), sup (1 — a)F _1 (a) > , 

« + 3 [a6(0,l/2] «e(l/2,l) J 

then r(l/n) = 0(n 2 / d {\og n) (d -^/ d ), whp. 

This mixing bound for GTG matches the best known mixing bound for RGG, and we believe 
that this equivalence is essentially correct. Intuitively, the an high weight nodes of a GTG G 
contain a spanning subgraph G' that is an RGG. The mixing time of G corresponds to the mixing 
time of G' . There are some extremely long edges in G\G' , but they seem to be too sparse to aid in 
mixing. At the same time, we find that the (1 — a) low weight nodes (with very short edges) do not 
slow down mixing. For technical reasons, we consider a weight decay of P[W > x] = 0{l/x d+u ). 
We conjecture that this equivalence continues to hold for ¥[W > x] = o(l/x l+l ^ d ). 



1 We will use the notation "with high probability" and denote it as whp, meaning with probability 1 — o(l) as n 
tends to infinity. 
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The rest of the paper is organized as follows. In Section [3l we derive upper and lower bounds 
on the maximal weight. Consequently, we find upper and lower bounds on the node degrees, and 
show that the number of edges |.E(Cr)| = O(ralogn). In SectionH]we construct a family of canonical 
paths for the GTG and then prove Theorem [2] in Section In Section [6j we reflect on our results 
and explain why we believe that our result holds for more slowly decaying weight distributions. 
Finally, in |Aj we exemplify our results with two different weight functions. Our first example is 
the exponential weight distribution f(w) = e~ w , with cumulative density function F(x) = 1 — e~ x . 
Our second example is the Pareto distribution with cumulative density function F(x) = 1 — x~ 7 , 
where x > 1 and 7 > d > 2. 



3 Node weights and node degrees in GTG 



In this section, we determine the upper and lower bounds on the maximal weight l^ mM in a 
geographical threshold graph G = (V(G),E(G)). Subsequently, we derive the upper and lower 
bounds on the degrees of the nodes in GTG near the connectivity threshold. Finally, we show that 
the number of edges |2£(G)| = O(nlogn) for these connected GTGs. 

We adopt the following notation for the remainder of the paper. We have a constant a G (0, 1) 
and we fix small constants e, v so that 



< e < v/2d. 



(5) 



Furthermore, we assume that there is a weight Wq such that if W > Wo then F[W > x] = 
0(l/x d+u ). For brevity, we will state this as "F[W > x] = 0{l/x d+u )" Finally, we use n(n) to 
denote an arbitrarily slowly increasing function of n, that is k(ji) = cj(1). 

The maximal weight W max satisfies P[W m ax < x] = F(x) n , since the weights are independently 
distributed. Consider a continuous weight distribution f(w) with cdf F(x). Our goal is to find two 
thresholds Wi,W 2 , such that P[W max < W x ] = o(l) and P[W max > W 2 ] = o(l). We can 'invert' 
F(x) using the quantile function F~ 1 (p) = inf{x G R + : p < F(x)}. Define 



Wi=F 



-1 



1 



n(n) 



n 



and W 2 = F- 1 [l 



1 



n n(n) 



We have P[W max < W{\ = (l- K (n)/n) n < exp(- K (n)) -> and P[W max < W 2 ] = (1- 
exp(— 1/k(ti) + 1/n K(n) 2 ) — > 1. In conclusion, the maximal weight VF max satisfies 

lim p[iy max e {W X ,W 2 )\ = 1. 



l/nn(n)) n > 



(6) 



See |A] for concrete examples of the calculation of the bounds on the maximal weight. 

Let us determine the upper and lower bounds for the node degrees, keeping in mind that the 
weight distribution has finite mean and variance. We consider the GTG around the connectivity 
regime, as described in Theorem [TJ The next result generalizes Lemma 3 in [16], which shows that 
all degrees of RGGs near the connectivity threshold are O(logn) whp. 



Lemma 3. Let G be a connected GTG with threshold function 6 n = cn/logn. 
v G V(G) satisfy deg(v) G Igtg where 



Igtg 



c\ logn, c 2 F 1 



1 



n n(n) 



logn 



Whp, the nodes 



(7) 



5 



for any function n(ri) = w(l) and for constants 



ci = I 1 



2c \ 2T d 
and C2 = . 



The minimum degree of a GTG is fi(logn). The maximum degree depends upon the decay 
rate of the weight distribution: a slower decay rate results in larger maximum degree. In El we 
calculate Igtg for weight distributions with exponential decay and polynomial decay. 

Proof. Using equation (j4j), we apply the Chernoff bound on the degree deg(v\w) of a node v with 
a given weight w: 

r deg(v\w) < (1 - <5)E[deg(u|w)]l < exp(-E[deg(v\w)]5 2 /2). 



Let C3 = 1 — yj < 1- By having E[deg(u|io)] = (n — l)p(tu) 
(5 = 1 — 03/(1 + w/fi), it follows 

deg(v|u>) < ci log n 



<exp( -^.(l + w/n) ( 1 



T d (n-1) 



C3 



w + fj,) log n and choosing 



1 + 



1 I logra 

n , 



(8) 



Next, we specify conditions such that equation is o(l/n) for all w > and sufficiently large 
n. For the sake of simplicity, let us denote x = l + w/[A>l, and consider the function 

c V x / 

The minimum of is attained at x = C3. Moreover, ^>(x) is strictly decreasing on (0,03) and 
strictly increasing on (03, +00). Because C3 < 1 and (f)(1) = — C3) 2 = 2, we know that <p(x) > 2 

for x > 1. That is, equation ([8]) is o(l/n), for n > 3. Thus, the degree distribution satisfies 



deg(-u) < c\ log n = J f(w)F deg(v \w) < c\ log n 
= o^n^ 1 / f(w)dw S j = o(l/n). 



dw 



The union bound gives the lower bound on degree of the nodes in the graph. 

We now obtain the upper bound. Equation ([6]) ensures F~ l (1 — «(n)/n) < W max < F^ 1 (1 — l/nn(n)) . 
Moreover, by the continuity of F~ 1 (x), for any e > 0, there is sufficiently large n = n(e), such that 
the upper and lower bounds on V^niax are arbitrarily close < Wi — W\ < e. 

The degree of the node with maximal weight satisfies the binomial distribution Bin(n— 1, (T^/^ n )(M / max + 
fj,)), which is concentrated around its mean (T^/c)(l — l/n)(T4 / max + /i)logn. Finally, the union 
bound gives the upper bound on the degrees. □ 

We now partition the interval Igtg of equation ([7]). We use this partition to calculate the 
number of edges | and again in Section [5] to bound the mixing time. 

Define h(x) = (1 — -F(x)) -1 , or equivalently, F(x) = 1 — l/h(x). This is also equivalent to 



h-\y)=F-\l-l/y). 
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By assumption we have h(x) = £l(x d+u ), so that h 1 (x) = O (x 1 /^ - ^) . Our first interval B 
contains the low weight nodes: 

B = {v G V{G) | w v < F~ x (l - a)} = {v G V{G) \ F{w v ) < 1 - a}. (9) 



Next, we partition the an nodes with weights in [F x (l — a), W max ]. By equation ©, VFmax 
F- 1 (1 - l/n/<n)) = r^n^n)). Let 

«0 = ^max 

a fc = maxji^ 1 (l - a^j^A ^"H 1 - «)} , for k > 1. 

The a/j are only defined until we reach — a). Call this final index M. Our partition consists 

of the subintervals of the form (a k ,a k -i\ for 1 < k < M. Note that the indexing of our endpoints 
is the reverse of the standard convention. 

Lemma 4. The final index satisfies M = o(log n) . 

Proof. We have a = W max = h^in^n)) = O ((nsfn)) 1 /^")). Let = (1 + e)/(d + v) < 1. By 
induction, for < k < M, 

a k = F- 1 (l - a-^) = = O (wQ = O ((«( n fN) . 

By definition, om = .F _1 (l — a). Let iV be the smallest integer such that {nK{n))^ N /( d+z/ ) < 
F -1 (l — a). The right hand side is a constant, so M = O(N) = 0(log log(nK(n))) = o(logn). □ 

For 1 < k < M, let 

A k = {v£ V(G) | w v G (a fe , a fc _i]} . (10) 
The degree of v G A k with weight u>„ is 

deg(f) = 0((1 + u>„)logn) = 0(a k -\ logra). (11) 
We now show that the sizes of the A k are all concentrated around their means. 
Lemma 5. Whp, for all 1 < k < M, we simultaneously have \A k \ = O {n/a^^j . 

Proof. We have F[W > a k ] = l/a\t\. Therefore E[|4i|] = n/aj +e and E[\A k \) = n/a l k t\-n/a\t e 2 = 
(n/alt{) for 2 < k < M. 

As for the concentration of these values, we consider the most delicate case of A\. We have 
F[v G A\] = P[w v > ai] = a ^ 1+t ^ = W m ix +t ^ ■= ?i- We can consider as being generated from 
Bin(n,gi). The Chebyshev inequality gives 

_ „ 91 | > „ ?l] < J_ = o f ^"»" +<l/ ' d+ ^ = 0( „-./ 2) . 

L J ngi \ n I 

Similarly, the probability that each remaining \ A k \ is more than twice its mean is 0{n~ 1 / 2 ). There 
are M = o(logra) such sets, so taking a union bound shows that all of them are concentrated 
whp. □ 

See [A] for concrete examples of the partition of the nodes according to weight. 
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Lemma 6. £^0 ^ = @U)- 

Proof. We use the general d'Alembert's convergence criterion: The sum of positive terms ^2k=o c k is 
convergent if there exist a positive integer N and 77 > such that k > N guarantees c k /c k _\ < 1— 77. 

Let 6 = F' 1 ^ ~ a) and b k = /i(6 A ._ 1 ) 1 /(i+^) for 1 < k < M. We have b k < a M -k-i for 
< k < M. It follows that 

y-< — + yl<i + yl. 

k=0 k M k=0 k k=0 k 

and 

1 /K _ fh-iY _ ( h-i 



V h ) \h(b k _iyni+<) 



< 1 - 77 



for large k since h{x) 3> x d+u S> x 1+e . Therefore this sum converges to a constant (independent of 
n). □ 

Lemma 7. If G is a GTG at the connectivity threshold then the number of edges \E(G)\ = 
0(nlogn) whp. 

Proof. Lemma [3] guarantees that all nodes have degree fi(logn) whp, so the number of edges 
= Sl(nlogn). As for the upper bound, whp 

M M 

2\E{G)\ = £ deg( W ) + £ £ deg(v) = o((|5| + £ LA fc |a fe _i) logn 

veB k=lv&Ak k=l 

M 

= Of an + ^2~[+^a k ^i) lognj = O(nlogn), 



k=l a k-l 



by Lemma [6l □ 



4 Canonical paths for GTG 

We employ canonical paths (as introduced in [25]) to calculate our bound on the mixing time. 
In this section, we construct the canonical paths for G, a connected GTG with threshold function 
9 n = cnj log n. For every ordered pair of nodes u,v E V(G) we choose a canonical path j uv between 
them. We define 

p = max - — S~] 7r(u)7r(7;)|7 w | 



max 



deg(u)deg(7;)|7™| , (12) 



e={x,y}GE(G) 2\E(G)\ 

juv 3>c 

where \^ uv \ is the length of the canonical path from u to v. As per [25] Proposition 12.1, the mixing 
time from node i satisfies 

n(5) ^(logvr^+logr 1 ), (13) 
where 6 > 0. We will set 5 = l/n, so that Tj(l/n) < p (log7r(7) _1 + logn) . 
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Let G be a GTG at the connectivity threshold 9 n = cn/logn for c < sup Qe ( ,i) a)/4, as 

per Theorem [TJ The proof of this result in [TTJ establishes the following two facts. First, a constant 
fraction of the nodes an have weights greater than — a). We let H(G) = {i £ V(G) \ Wi > 

— a)} denote this set of high weight nodes and let L(G) = {i £ V(G) \ Wi < — a)} 

denote the complementary set of low weight nodes. Second, each high weight node is connected to 
every node within the critical radius 

log(an) \ 1/d 
anT d J 

In other words, the induced subgraph on H(G) contains a subgraph G' that is a connected RGG. 
Every vertex in L(G) is adjacent to nodes in H(G), so G is also connected. 

We use the connected RGG subgraph G' to construct our canonical paths. Our construction for 
the canonical paths is similar to the one used in [16] to bound the mixing time of RGG. Compared 
to that result, our proof addresses a novel technical challenge: all of the degrees of a RGG are 
O(logn), while the degrees of a GTG are heterogeneous. This leads to two challenges. First, when 
u, v E L(G), the intermediate nodes on the canonical path 7 U1) must all be in H(G). Therefore, 
we must choose these paths so that the intermediate nodes are evenly distributed among the high 
weight nodes. Second, many high weight nodes have degrees that are cj(logn). We must ensure 
that their contribution to the sum in equation (|12p does not lead to an increase in the mixing time, 
compared to a RGG. 

We now describe the geometric scaffold for our canonical paths, as in [16]. Partition the unit 
cube into a toric grid of [k] d small cubes, where k is specified below. A set of canonical paths for 
the grid will act as the framework for our canonical paths for the GTG. Intuitively, we increment 
the entries (mod k) in succession. So first we increase a\ until we achieve b\, then do the same for 
the second entry, and so on. The canonical path from (oj, a%, ■ • • , ad) to (61, b 2 , ... , b d ) is 

(01, 02, • • • , a d ), (01 + 1, a 2 , ■ ■ ■ , a d ), • • • , (h - 1, a 2 , ■ ■ ■ , ad), (bi,a 2 , ... , a d ), 
(bi,a 2 + 1,. . .,a d ), (h,b 2 , ... ,b d -l), (&i,6 2 , • ■ • ,b d ). 

Each path has length at most dk. Note that we always increment the index by +1 (even if there is 
a shorter path). 

While there are k 2d canonical paths, each edge appears in no more than k d+l paths. In- 
deed, any path that includes the edge from (ii, . . . i t , . . . i d ) to . . . it + 1, . . . i d ) must start at 
{h,... l t ,jt+i,- ■ ■ ,jd) and end at (jx, . . . j t - 1, l' t , ■ ■ ■ , l' d ) for some l\,...l t and l' t , . . . , l' d . This results 
in t choices for 1%, . . . , If and d — t + 1 choices for l' t , . . . l' d . 

Before constructing the canonical paths for our GTG, we must prove two lemmas. Tile [0, l] d 
into cubes Si with side length l/k = (c' log n/n) l / d = 0(r conn ), where we state the conditions on 
constant c' later. We have a [k] d grid, whose cubes each have volume c'logn/n. Let H(Si) and 
L(Si) denote the high weight and low weight nodes in Si, respectively. 

Lemma 8. There exist constants /3q, fii > such that whp every cube S satisfies /?o logn < ^(.S)! < 
Pi log n and (3q log n < \H(S)\ < (3i log n. 

Proof. This proof is similar to the proof of Theorem 3 in [11] . and will use the same notation. Let 
Bi = \H(Si)\ and Ri = \L(Si)\. In expectation, there are E[2?j] = ac'logn high weight nodes within 
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Si. Using the lower and upper tail Chernoff bounds [151 . it follows 



% 1(1- <h }EM 1 (<xp ( -jE[Bi] ) . 



' [Bi > (1 + 5 2 )E[B l ]] < exp 



61 MBi 



2(l + <5 2 ) 



Fixing 62 G [0,1], we take ft = ft/ VI + <^2 , and thus ft G [0,1]. The number of high weighted 
nodes Bi within the cube Si satisfies 



' [Bi G (1 — ft, 1 + ft) E[Bi]] > l-2exp 



51 



2(1 + ft) 



-E\Bi 



1 - 2n 



-ac'6%/(l+8 2 ) 



By the union bound it follows: 

P g (1 - ft, 1 + ft) ![£;]} 



> 1 - 2n" 



c'«52/(l+5 2 )^ 



n/(c' log n) 



Taking the limit of the last expression as n — > +00, for d > (l+ft)/a;5|, we obtain the concentration 
on Bi, for each cube Sf. 



lim I 

n— >oo 



f|{Si G (1 - ft, 1 + ft)ac' log n} 



The concentration on the number of low weight nodes Ri, within each cube Si, follows analogously 
to the previous analysis. In expectation, there are E[i?j] = (1 — a)d log n low weight nodes within 
Si . Hence by the Chernoff tail bounds and the union bound we have 



lim P 

n— >oo 



f|{i?i G (1 - ft, l + 5 2 ) (!-«)</ log n} 



Finally, we can guarantee the concentration of both the high and low weight nodes by taking 
c '> max {^ ; _^},thatis, 



d> 



^, if a G (0,1/2) 
i+s 2 



if a G [1/2,1). 



(16) 



So the lemma holds with /3o = (l — ft/Vl + ft) min{a, 1 — a} and /ft = (1 + ft) max{a, 1 — a}. □ 
Two cubes are adjacent if they share a (d — l)-dimensional boundary. 

Lemma 9. Let Si, Sj be adjacent cubes. The number of edges between H(Si) and H(Sj) is f2(log 2 n) 
provided that c < (d + 3) -d / 2 min{sup a , g ( 1 / 2i i)(l — x)F~ 1 (x), sup^g^ 1/21 xF^ 1 (x)}. 



See [A] for explicit calculation of the constant c for two example weight distributions. 
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Proof. Consider any high weight node u G H(Si) and any high weight node v G H(Sj), with the 
weights w u and w v , respectively. The distance r uv between u and v is at most \/d + 3(c' log njn) x l d . 
Indeed, the furthest points in two adjacent d-dimensional unit cubes are at distance y/d~+H by the 
(i-dimensional Pythagorean theorem. 

Consider the connectivity relation 

Wu + Wv F^il -a)+ F'ijl -a) _ 2F~ 1 (1 - a) n 
ri v ~ (d + 3) d / 2 c' log n/n ~~ (d + 3) d / 2 c' logn" 

High weight nodes u,v are connected with probability one if (w u + w v )/r d v > 9 n = cn/logn, which 
is guaranteed if 2i ?_1 (l — a)/((d + 3) d ' 2 d) > c. Using equation (JT6J) , we require 

2 b\ f qf-^l-q), if a G (0,1/2) 

° " (d + 3) d / 2 1 + & I (1 - ajF-^l - a), if a G [1/2, 1). 

Since 52 G (0, 1) is arbitrary, with sup^g^ ^ 5 2 /(l + 62) = 1/2, the conditions combine to give 
c < (d + 3) _c!/2 min{ sup a.F~ l (l-a), sup (1 - a)F _1 (l - a)) 

06(0,1/2) ae[l/2,l) J 

= (d + 3) _d/2 min{ sup (l-x)F _1 (x), sup xF _1 (x)). 

se(l/2,l) ae(0,l/2] J 

□ 

We now employ a randomized procedure for choosing canonical paths. This procedure guar- 
antees that no edge appears in more than Q(k d+1 ) = OfV^" 1 ) = 0((n/ dog n)^ d+l ^ d ) canonical 
paths. Let us denote cubes as S(a±, . . . , a^) for 1 < a» < k, where (ai, ... ,0,4) is the location 
in the d-dimensional grid. By Lemma EJ we have both /?ologn < \L(S(ai, . . . ,ad))\ < /3ilogn 
and /3ologn < |i7(5(ai, . . . , a^))| < /3ilogn for some constants /3o, /?i > 0. For each cube 
S = S(ai, ■ ■ ■ , ad), evenly partition L(S) into sets Li(S), for 1 < i < /3ologn. Each set has 
size at most [73i/A)l whp. Assign each set of low weight nodes to a distinct high weight node 
hi G H(S), for 1 < i < /3ologn. The high weight node hi is the high weight representative for the 
nodes in Li(S). The key outcome of this assignment is that each h G H(S) represents a constant 
number of low weight nodes whp. 

Consider any ordered pair of nodes (x,y). Let x G S(a\, . . . , oy) and y G S(bi, . . . , bd)- We 
choose a canonical path from x to y as follows. We use the toric grid to identify the sequence of cubes 
in the canonical path. Taking equation (fT5|) as our framework, we consider cubes S(ai, 02, • • • , ad), 
S(a 1 + l,a 2 ,. ■ ■ ,a d ), S(h, a 2 , ■ ■ ■ , a d ), S(bi,a 2 + l, ■ ■ ■ ,a d ), ■■■ , S(b 1 ,b 2 , ■ ■ ■ ,b d ). For brevity, call 
these cubes Sq, S±, . . . , St- If x is a low weight node, then set xq to be the high weight representative 
for x. If x is a high weight node, set xo = x. For 1 < i < t, choose Xj to be a random high weight 
node in Si- Set xt = y if y is a high weight node, otherwise using the high weight representative 
for y. Our canonical (x, y) path is x, xq, x±, . . . , xt, y. If xq = x or xt = y, we remove the repeated 
node from the path. We have 




and furthermore, we can bound how often each edge appears in a canonical path. 
Lemma 10. Every edge in G appears in at most 0((n/ log n)( d+1 ^' d ) canonical paths whp. 
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Proof. Let Z xy denote the number of times the edge xy is chosen. If both x and y are low weight 
nodes then Z xy = (we always move to a high weight node from a low weight node). When x is 
low weight and y is high weight, then the edge xy can only be used when x and y are in the same 
cube. In this case, the edge xy can only be chosen when the canonical path has x as one of its 
cndpoints. Therefore Z xy < 2(n — 1). 

Now suppose that both x and y are high weight nodes. The edge xy is only used if x, y are in the 
same cube or in adjacent cubes. First, we consider high weight x, y in the same cube S. The edge xy 
will be used only by paths between Lj(S') U {x} and Lj(S) U {y}. If one or both is not a high weight 
representative, then this edge will be used even fewer times. Therefore Z xy < (1+ \f5\/ (5q\) 2 = 0{\) 
whp. 

Consider x G S and y G T in adjacent cubes 5, T. Let 7^ = {7 U1 , | xy G 7™}. We consider 
four cases, according to the locations of u and v. If u G S and v £ T then just as in the previous 
case, whp there are at most (1 + |7?i/A)l) 2 = 0(1) such paths. 

The total number of paths j uv with u G S and v ^ T is 0(n). Indeed, x must be the high 
weight representative of u, which gives 1 + \Pi//3o] = O(l) choices for u, and there are 0(n) end 
nodes for paths that start in S\. Similarly, the total number of paths with u ^ S and v G T is 
0{n). 

The remaining case is when u G" 5 and v ^ S. Let denote the number of times the edge xy 
is chosen as the random edge from S to T for some canonical path. Let 

V(S, T) = {^f uv I 7 Mt) traverses from S to T and u ^ S andv ^ T}. 

Then whp 

2 1 _ /" AV„-<i-i 



E[Zg = £ P(xy G 7fW ) < r-^Vi log n) 2 — — = 

~/uv£'P(S,T) 

Indeed, the number of canonical paths for the toric grid which pass from S to T is CK^nn 1 )- 
Each edge between cubes corresponds to at most logn) 2 toric paths in the GTG. Since both x 
and y are internal nodes of these paths, they were each chosen uniformly and independently with 
probability at most 1 / (j3o log n) . 

Using the Chernoff bound for this binomial distribution, 

2 

P [\Z' xy - E[Z' xy ]\ > r)E[Z xy \] < 2exp ( - |~E[^]) 

- M4<f> 2 (^> ,d+,,/d )^">. 

The union bound now gives 

F[A Xty (\Z' xy -E[Z' xy }\> v E[Z' xy })] 
< £P [\Z' xy - E[Z' xy }\ > V E[Z' xy }] < n 2 o(l/e") 0. 



<i-l 
1 1 j ■ 1 r 



Therefore whp, every edge between high weight nodes in adjacent cubes is used by (1 ± r/)r l 
e((n/logn)( d + 1 )/ c( ) canonical paths. In this case, E[Z xy ] = 6((n/ log n)^ d+l ^ d ). □ 
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5 The mixing time for GTG 



In this section, we prove that the mixing time for a ci-dimensional GTG near the threshold for 
connectivity is 0(n 2 / d (log n) i - d - 2)/d ) when F[W > x] = 0(l/x d+u ) for v > 0. 

We have | i?(C) | = ©(nlogn) and max u ^, g y( G ) \j U v\ = 0((n/ log n) l / d ) by Lemma[7]and equa- 
tion (fT7|) . respectively. Substituting these values into equation fp~2|) yields 

p < O ( 1W ,, — -7T-- 7TTT | ( max deg(-u) deg(w) I . (18) 

W^V^logn) A* J \ e={x,y}£E(G) ^ &K ' &V ' ) V ' 

Fix an edge e = {x, y} between high weight nodes in adjacent cubes, and define cr(e) = ^ 7uij9e deg(it) deg 
Suppose that these high weight cubes differ in the kth coordinate. Specifically, the cubes contain- 
ing x and y are indexed by (ji, ■ ■ ■ , jk-i, jk, jk+i, ■ ■ ■ , jd) and (ji, . . . Jk-iJk + lJk+U---Jd), 
respectively. If e G 7™ then u must be in an initial cube indexed by (ii, . . . , ifc-i, ifc> Jfe+ii ■ ■ ■ ,jd), 
while v must be in a target cube indexed by (ji, . . . , Jk-l Aki^k+Xi ■ ■ ■ ^d)- Let U\ C [0, l] d corre- 
spond to the set of initial cubes, and let U2 C [0, l] d correspond to the set of target cubes. The 
volumes of these subsets are vol(t7i) = [nj log n) k / d /(n/ log n) = (\ogn/n)^ d ~ k ^ d and woKJJ^) = 
[nj log n)( d ~ k+l ^ / d / (n / log n) = (log n/n) <yk ~ 1 ^ d . Let V\ (respectively V2) be the set of nodes in U\ 
(respectively U2). 

In order to bound er(e) we use the partition of the vertices into the sets B and A%, A2, . . . Am 
as defined in equations ([9]) and (jlOp . We consider canonical paths containing e that run between 
every possible pair of sets in this partition. The following three technical lemmas require the weight 
decay of ¥[W > x] = 0{l/x d+u ). 

Lemma 11. For anyT e {B, Ai,A 2 , . . .A M }, let Zi = \TC\Vi\ fori = 1,2. ThenK[Zi] = \T\vol(Ui) 
and there exists a constant C > ( depending on T and i ) such that 

P [\Zi - E[Zi]\ > E[Zi\] = O (ex 9 {-Cn v/2d{ - d+v) ) 
In other words, the probability that Zi ^ 0(E[Zj]) decays at a super-polynomial rate. 

Proof. We prove this concentration for Z = Am H V\, which is the most delicate case. The other 
cases follow similarly. First, by the independence of weights and location, we have \Am H V\\ 
Bin(|AAf|, vol(Lq)). Furthermore, by Lemma \Am\ = ©(n/aj^). Therefore Am H V\ 
BmiC'n/a 1 ^, (logn/n)^-^^) for some C > 0. The expected value is 

C n /log n\ C' n flogn^^ d ~ k '' d 



m = -t+t -z- > 



> 



O^M-l V n J WAtx V II 

C'n (logra) 1 ^ 



(nK(n))( 1 + £ )/( rf +^) n^- 1 )/^ 

> C ' n yd-(l+e)/(S+u) = n ^ n u/2d(d+u^ (lg) 

since e < u/2d by equation ([5]) and «(n) grows arbitrarily slowly. Using the binomial Chernoff 
bound (cf [Mj, Corollary 2.3), we have 

¥[\Z-E[Z}\ > E[Z\] < 2exp(-^K[Z]^ = O (eicp{-Cn v/2d{d+u) ) 

for some constant C > 0. □ 
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For T\, T 2 G {B, A±,A 2 , ■ ■ ■ Am}, let A(Ti, T 2 ) denote the number of canonical paths from Tin Vi 
to T 2 Pi V 2 that use edge e. 

Lemma 12. ForTi,T 2 G {5, Aj, A 2 , . . . Am}, we /iave 

E[ A ( r,r 2 )i = e(^M^M) 

V log n J 

and there exists a constant K ( depending on Ti,T 2 ) such that 

P[A(Ti,T 2 ) -E[A(Ti,r 2 )]| > E[A(Ti,r 2 )]] = O (exp(-Kn^ d+ ^ . 

Proof. The proof is similar to the previous one. For % = 1,2, let = T,nVi, so that \Z{\ = @(E[Zj]) 
by Lemma 1111 The distribution for A(Ti,T 2 ) is Bin(|Zi||Z 2 |, K'/log 2 n) for some K' > 0, and 
\Zi\ = 9(E[Z<]) for i = 1,2 by LemmaEl 

Furthermore, \Z { \ = ^n"/ 2 ^ 4 "^) by equation flS). Therefore E[A(T X , T 2 )\ = n(n u / d( - d+u ^ / log 2 n) = 
UJ ^ n u/2d(d+u)y rp^g kj nom j a i Chernoff bound gives 

P[|A(ri,r 2 )-E[A(Ti,r 2 )]|> E[A(ri,r 2 )]] < 2exp(-~E[A(ri,T 2 )]) 

= O (expi-Kn^^fj 
for some constant K > 0. □ 



The previous two lemmas show that these quantities are tightly concentrated around their 
means. A union bound shows that whp these quantities are concentrated for all candidate edges 
simultaneously. Indeed, there are 0(d ■ (n/ log n) 1 ^) choices for adjacent cubes used in canonical 
paths, with 0(log 2 n) edges running between each pair. Recall that M = o(logn) by Lemma [H 
For a fixed choice of adjacent cubes, there are 2(M + 1) = o(logra) choices for T in Lemma [TTJ 
Our union bound for these event involves 0((n/ log n) l l d log n) terms, each decaying at a super- 
polynomial rate. Therefore, all these events are concentrated whp. As for Lemma [I2l there are 
(M+ l) 2 = o(log 2 n) choices for (Ti, T 2 ) for each pair of adjacent cubes. This union bound is taken 
over all relevant edges between high adjacent cubes. The number of addends in this union bound is 
0((n/ log n) 1 /^ (log 2 n)(log 2 n)). Again, the concentrations from Lemma [T2l are super-polynomial, 
so the union bound shows that all these quantities are concentrated simultaneously whp. 

Lemma 13. Our choice of canonical paths gives p = 0((logn) 2 / d ). 

Proof. We bound a(e) = 2~ 9e deg('u) deg(u ). First consider the contributions from canonical 
paths from B f] V\ to B n V 2 . Let a(B,B) denote the contribution of paths between low weight 
nodes to cx(e). 

Recall that if v G B then deg(f) = 0(logn). By Lemmas 1121 and 111} we have whp 

E[a(B,B)} = o(^^--(log 2 n)) = O (\B n V x \ ■ \B n V 2 \) 
V log n / 

_/ {\ogn\ ( ' d ~ k ^ d flogn^ 
O n • n 



n I \ n 



o{n^l d {\ogni d -^l d ) 
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We consider the contribution of paths between high weight nodes by using the partition A±, A2, ■ ■ ■ 
specified in equation (fTO]) . Let a(A r ,A s ) denote the contribution to a(e) for paths from A r to A s . 
Recall that, as per equation (fTTI) . if v 6 A^ then deg(f) = 0(ak-i logn). Arguing similarly to 
the calculation above, the contribution to cr(e) from paths between nodes in A r and A s , where 
< r, s < M, is whp 

a(A r ,A s ) = O ( A( i yLr ; Ag) • a r a s log 2 n 
V log n 

( n /\ogn\( d ~ k ^ d n (\ogn\^ k ~ l ^ d 



/ n (rf+l)/rf(l ogn )('i-l)/'i' 



Next we consider the paths between low weight and high weight nodes: whp 

log 2 n 



a(B,A s ) = o(^V^-a s log 2 n 



'logn\^ rf k ^ d n /logn\^ 



o 



and similarly, a(A r ,B) = O (n^ d+1 ^ d (logn)^ d l '' d /a^). Putting these estimates together, whp 
every edge e between cubes satisfies 

M M M 

(7(e) = v(B, + Aj) + a(Aj,B)) + J^Y, °( A i> A *) 

3=1 3=1 k=l 

(M MM 
n^(logn)( d - 1 )/-(l + 2X;i + EE^ 
i=i « j=i fc=i j fc 



(M M \ 

n (d+D/<i (logn) (d-l)/^ 1 + 2 ^i_ + ^J_^ ) 
t=l * j=l i J 



= 0[n^l d (\ogn)^/ d 
where the last equality follows from Lemma [6j 



Finally, by equation (|18p we have whp p = O ( (ra/logn) 2 ^ ) . □ 



Proof of Theorem[2l Equation (fT3|) gives r x (<5) < p (log 7r(x) 1 + log 5 x ) . The previous lemma 
ensures that p = 0((n/ log n) 2 l d ). Meanwhile, we have ir(x) = deg(x)/2\E(G)\ = $7 (logn/ (n logn)) = 
ri(l/n) by Lemma [3] and Lemma [71 In summary, r x (,J) = 0(n 2 / d (logn)( d ~ 2 )/ d ) for 5 = l/n. □ 
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6 Conclusion 



We have shown that if the weight distribution of a d-dimensional GTG satisfies ¥[W > x] = 
0{l/x d+v ), then its mixing time is 0(n 2 /d(log n)^ d 2 ">l d ). This matches the known bounds for 
RGG. Our proof uses a spanning subgraph among the an high weight nodes to create a scaffold 
for canonical paths. In constructing these paths, our proof treats all high weight nodes identically, 
ignoring the particularly large reach of the highest weight nodes. We did try to take advantage of 
these hub nodes, but found that they were to sparse to leverage for canonical paths. 

One might wonder whether this is a shortcoming of the method of canonical paths, rather than 
a reflection on the characteristics of GTG. However, initial investigations using conductance to 
bound mixing (as in [4J) suggests the same conclusion. Of course, using conductance introduces its 
own set of technical challenges due to the heterogeneous nature of the degrees. 

For technical reasons, we assumed that the weight distribution decayed as ¥[W > x] = \/x d+v . 
We believe that the equivalence of mixing for GTG and RGG extends to weight distributions with 
slower decay. In particular, we conjecture that these mixing time of a GTG with weight decay 
¥[W > x] = 0(1/ x 7 ) matches that of RGG when 7 > 1 + 1/d, and that GTG mixes faster when 
1< 7 < 1 + l/d. 
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A Characteristics of GTG for example weight distributions 

We describe the relevant characteristics of GTGs for two different weight distributions: exponential 
decay and polynomial decay. 

A.l Exponential Weight Distribution 

Our first example is the exponential weight distribution f(w) = e~ w with cumulative density 
function F(x) = 1 — e~ x . Inverting the cdf gives F~ 1 (x) = — log(l — x). 

We first discuss the weights and degrees of the nodes in the GTG, as described in Section [3j 
As per equation (|6|), the maximum weight satisfies 

W max G [log n - K(n), log n + n(n)} 

whp. Lemma [3] guarantees that whp all the node degrees are in the interval 

Igtg= [O(logn),0(log 2 n)] . 

Next, we partition the interval Iqtg- By equation ([9]), the cutoff for low weight nodes is 
F _1 (l — a) = log a -1 . We partition the high weight nodes into disjoint subsets as specified in 
equation (TlOl) . Note that F^ 1 (l — x~( 1+e )) = log(x 1+t ) = O(logx). Therefore, ignoring leading 
constants, the sequence of endpoints (in descending order) is 

( log n, log log n, log log log n, . . . , log^ M ^ nj 

where M < log* n, where the iterative logarithm function log* n is the number of iterations of the 
log function required to obtain a result less than 1. 

Finally, we calculate upper bound on the constant c required by Lemma [9] in Section [4l For the 
exponential distribution, taking c < ((d + 3) d//2 e) _1 is sufficient. 
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A. 2 Pareto Weight Distribution 



We now give a parallel characterization of our second example: a Pareto distribution with cumula- 
tive density function F(x) = 1— x~ 7 where 7 > d > 2. Inverting this cdf gives F~ 1 (x) = (1— x)^ 1 ^ . 
We consider our results concerning node weights and node degrees. We have 



K\n) 



1/7 



(n K (n)) 1/7 



whp by equation ([6]). By Lemma El whp all the node degrees are in 

Igtg = [6(log n^ean^n)) 1 ^ logn)]. 

We separate the low weight nodes from the high weight nodes using the weight cutoff — a) = 

afVT, Next, we partition the high weight nodes as per equation (jlOp . Note that F^ 1 (l — x~( 1+e )) = 
x( 1 + e )/')' - it follows that our sequence of endpoints is 

((n K(n)) 1 ^, (n n(n))^, (n n{n)f 2 ^, . . . , (n K {n)f U ^) 

where (3 = (1 + e)/7 < 1. Here M is the smallest integer such that (nn(n))^ ^ < .F _1 (l — a) = 
a -1 / 7 . The latter requirement is equivalent to the condition M > (loglog(n«(n))— loglog(a _1 ))/log(/3 -1 ) 

Finally, we calculate an upper bound on the constant c required by Lemma [9] in Section 01 Since 
F^ 1 (x) = (1 - x)- 1 ^, by Lemma[9]it follows c < {{d + 3) d / a 2(7-i)/7)-i. 
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